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Abstract 

The isospin asymmetries of i? — ?■ K^fi + fi~ decays and the partial branching frac- 
tions of B° — > K°[i + {i~ and B + — > K* + fi + [i~ are measured as a function of the 
di-muon mass squared q 2 using an integrated luminosity of 1.0 fb _1 collected with 
the LHCb detector. The B — > K[i + [i~ isospin asymmetry integrated over q 2 is neg- 

i—t ative, deviating from zero with over 4 a significance. The B — > K* /i + fi~ decay mea- 

surements are consistent with the Standard Model prediction of negligible isospin 

"k> asymmetry. The observation of the decay B° — > K®[i + fi~ is reported with 5.7 a 

significance. Assuming that the branching fraction of B° — > K°[i + fi~ is twice that 
of B° ->■ K®fi + fj,~, the branching fractions of B° -)■ K° /i + fi~ and B ->■ K* + fi + fi~ 
are found to be (0.31^;^) x 10~ 6 and (1.16 ± 0.19) x 10~ 6 , respectively. 
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i 1 Introduction 



2 The flavour-changing neutral current decays B— y K^fi + fi~ are forbidden at tree level in 

3 the Standard Model (SM). Such transitions must proceed via loop or box diagrams and 

4 are powerful probes of physics beyond the SM. Predictions for the branching fractions 

5 of these decays suffer from relatively large uncertainties due to form factor estimates, 
e Theoretically clean observables can be constructed from ratios or asymmetries where the 
7 leading form factor uncertainties cancel. The CP averaged isospin asymmetry (A{) is such 
s an observable. It is defined as 



9 where T(B — > /) and BiB — > f) are the partial width and branching fraction of the B — > f 

10 decay and tq / t + is the ratio of the lifetimes of the B° and B + mesons Q For B — > K* [iT , 

11 the SM prediction for A\ is around — 1% in the di-muon mass squared (q 2 ) region below 

12 the J/ip resonance, apart from the very low q 2 region where it rises to 0(10%) as q 2 

13 approaches zero [1]. There is no precise prediction for A\ in the B— > Kfi + fi~ case, but it 
u is also expected to be close to zero. The small isospin asymmetry predicted in the SM is 

15 due to initial state radiation of the spectator quark, which is different between the neutral 

16 and charged decays. Previously, A\ has been measured to be significantly below zero in 

17 the q 2 region below the J/ip resonance [21 E] • In particular, the combined B — > Kfi + fi~ and 
is B— > K*fi + fi~ isospin asymmetries measured by the BaBar experiment were 3.9 a below 

19 zero. For B—t K*y + fi~ , A\ is expected to be consistent with the B — > K* 0/ y measurement 

20 of 5 ± 3% jl] as q 2 approaches zero. No such constraint is present for B—t Kfi + yT . 

21 The isospin asymmetries are determined by measuring the differential branch- 

22 ing fractions of B+ ^ K+ fx~ , B° , B° -> (K*° K + n~)fi + /i~ and 

23 B + — 7- (K* + — > K®7T + )fi + {i~] the decays involving a K® or 7r° are not considered. 

24 The meson is reconstructed via the K® — > 7t + tt~ decay mode. The signal se- 

25 lections (Section [3]) are optimised to provide the lowest overall uncertainty on the 

26 isospin asymmetries; this leads to a very tight selection for the B + — > K + fi + and 

27 B° — > (K*° — > K + 7T~)/j + fi~ channels where signal yield is sacrificed to achieve overall uni- 

28 formity with the B°— > K®fi + [i~ and B + — > (K* + — > K®n + )fi + fi~ channels, respectively. 

29 In order to convert a signal yield into a branching fraction, the four signal channels are 

30 normalised to the corresponding B— > J/tpK^ channels (Section fsj. The relative normal- 

31 isation in each q 2 bin is performed by calculating the relative efficiency between the signal 

32 and normalisation channels using simulated events. The normalisation of B° — > K°fi + fi~ 

33 assumes that B(B° — >• K°fi + fi~) = 2B(B° — >• K®fi + fi~). Finally, Ai is determined by si- 

34 multaneously fitting the n + fi~ mass distributions for all signal channels. Confidence 

1 Charge conjugation is implied throughout this paper. 



Ai 



r(5°-+ K(*)°fj,+n-) + T(B+^ KW+fi+fi-) 
B{B -^ K^°fi + fi-) - ^B(B+^ K^+/j,+ix-) 



(1) 



B(B ^ K(*)°n+n~) + ^B(B+^ K^+n+fi-)' 
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35 intervals are estimated for A\ using a profile likelihood method (Section H). Systematic 

36 uncertainties are included in the fit using Gaussian constraints. 

37 2 Experimental setup 

38 The measurements described in this paper are performed with 1.0 fb™ 1 of pp collision data 

3 9 collected with the LHCb detector at the CERN LHC during 2011. The LHCb detector [5] 

40 is a single-arm forward spectrometer covering the pseudorapidity range 2 < 77 < 5, de- 

41 signed for the study of particles containing b or c quarks. The detector includes a high 

42 precision tracking system consisting of a silicon-strip vertex detector (VELO) surrounding 

43 the pp interaction region, a large-area silicon-strip detector (TT) located upstream of a 

44 dipole magnet with a bending power of about 4Tm, and three stations of silicon-strip 

45 detectors and straw drift-tubes placed downstream. The combined tracking system has a 

46 momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, and 

47 an impact parameter (IP) resolution of 20 |J.m for tracks with high transverse momentum. 

48 Charged hadrons are identified using two ring-imaging Cherenkov (RICH) detectors. Pho- 

49 ton, electron and hadron candidates are identified by a calorimeter system consisting of 

50 scintillating-pad and pre-shower detectors, an electromagnetic calorimeter and a hadronic 

51 calorimeter. Muons are identified by a muon system composed of alternating layers of 

52 iron and multiwire proportional chambers. 

53 The trigger consists of a hardware stage, based on information from the calorimeter 

54 and muon systems, followed by a software stage which applies a full event reconstruction. 

55 For this analysis, candidate events are first required to pass a hardware trigger which 

56 selects muons with a transverse momentum, pp > 1-48 GeV/c for one muon, and 0.56 and 

57 0.48 GeV/c for two muons. In the subsequent software trigger [£], at least one of the final 

58 state particles is required to have both p T > 0.8GeV/c and IP > 100 |xm with respect to 

59 all of the primary proton-proton interaction vertices in the event. Finally, the tracks of 

60 two or more of the final state particles are required to form a vertex which is significantly 

61 displaced from the primary vertices in the event. 

62 For the simulation, pp collisions are generated using Pythia 6.4 [7j with a specific 

63 LHCb configuration [8j. Decays of hadronic particles are described by EvtGen [9] in 

64 which final state radiation is generated using Photos [TU]. The EvtGen physics model 

65 used is based on Ref. [TT]. The interaction of the generated particles with the detector 
ee and its response are implemented using the Geant4 toolkit [12] as described in Ref. [13J. 

67 3 Event selection 

68 Candidates are reconstructed with an initial cut-based selection, which is designed to 

69 reduce the dataset to a manageable level. Channels involving a K® meson are referred 

70 to as K® channels whereas those with a K + meson are referred to as K + channels. 

71 Only events which are triggered independently of the K + candidate are accepted. 

72 Therefore, apart from a small contribution from candidates which are triggered by the 
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73 Kg meson, the K® and the K + channels are triggered in a similar way. The initial 

74 selection places requirements on the geometry, kinematics and particle identification 

75 (PID) information of the signal candidates. Kaons are identified using information from 

76 the RICH detectors, such as the difference in log-likelihood (DLL) between the kaon 

77 and pion hypothesis, DLLxtt- Kaon candidates are required to have DLL^ > 1, which 

78 has a kaon efficiency of ~ 85% and a pion efficiency of ~ 10%. Muons are identified 

79 using the amount of hits in the muon stations combined with information from the 
eo calorimeter and RICH systems. The muon PID efficiency is around 90%. Candidate 
ei are required to have a di-pion mass within 30 MeV/c 2 of the nominal K® mass 

82 and K* candidates are required to have an mass within 100 MeV/ c 2 of the nominal K* 

83 mass. At this stage, the K® channels are split into two categories depending on how 

84 the pions from the K® decay are reconstructed. For decays where both pions have hits 

85 inside the VELO and the downstream tracking detectors the K® candidates are classified 
se as long (L). If the daughter pions are reconstructed without VELO hits (but still with 
87 TT hits upstream of the magnet) they are classified as downstream (D) K® candidates, 
as Separate selections are applied to the L and D categories in order to maximise the 
89 sensitivity. The selection criteria described in the next paragraph refer to the K® channels. 

90 

91 After the initial selection, the L category has a much lower level of background than 

92 the D category. For this reason simple cut-based selections are applied to the former, 

93 whereas multivariate selections are employed for the latter. Both B° and B + L selections 

94 require the K® decay time to be greater than 3 ps, and for the IP x 2 t° be greater than 10 

95 when the IP of the K®, with respect to the PV, is forced to be zero. The B°— >■ K°fi + fi~ 

96 L selection requires that K® px > 1 GeV/ c and B > 2 GeV/ c. The K® mass window is 

97 also tightened to ±20 MeV/ c 2 . The B + — > K* + fi + ^ L selection requires that the pion from 

98 the K* + has an IP \ 2 > 30. Multi-variate selections are applied to the D categories using 

99 a boosted decision tree (BDT) [14J which uses geometrical and kinematic information of 

100 the B candidate and of its daughters. The most discriminating variables according to 

101 the B° and B + BDTs are the K® pi and the angle between the B momentum and its 

102 line of flight (from the primary vertex to the decay vertex). The BDTs are trained and 

103 tested on simulated events for the signal and data for the background. The simulated 

104 events have been corrected to match the data as described in Sect. HI All the variables 

105 used in the BDTs are well described in the simulation after correction. The background 
loe sample used is 25% of B candidates which have \m K ( i *) IJt + „- — m B \ > 60 MeV/ c 2 , where 
107 is obtained from fits to the appropriate B— > J/^K^ normalisation channel. These 
los data are excluded from the analysis. The selection based on the BDT output maximises 
109 the metric S/\/S~+B, where S and B are the expected signal and background yields, 
no respectively. 

in The K + channels have, as far as possible, the same selection criteria as used to select 

112 the K^ channels. The cut-based selections applied to the L categories have the K® specific 

in variables (e.g. K® decay time) removed and the remaining requirements are applied to 

H4 the K + channels. The BDTs trained on the D categories contain variables which can be 

us applied to both K® and K + candidates and the BDTs trained on the K® channels are 
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116 simply applied to the corresponding K + channels. The K + channels are therefore also 

117 split into two different categories, one of which has the L selection applied, while the other 
us one has the D selection applied. The overlap of events between these categories induces 
n9 a correlation between the L and D categories for the K + channels. This correlation is 

120 accounted for in the fit to A\. 

121 The final selection reduces the combinatorial background remaining after the initial 

122 selection by a factor of 5-20, while retaining 60-90% of the signal, depending on the cat- 

123 egory and decay mode. It is ineffective at reducing background from fully reconstructed 

124 B decays, where one or more final state particles have been misidentified. Additional 

125 selection criteria are therefore applied. For the channels, the A — > pir~ decay can be 

126 mistaken for a K® — > 7i + tt~ decay if the proton is misidentified as a pion. If one of the 

127 pion daughters from the K® candidate has a DLL p7r > 10, the proton mass hypothesis is 

128 assigned to it. For the L(D) categories, if the pjr~ mass is within 10(15) MeV/c 2 of the 

129 nominal A mass the candidate is rejected. This selection eliminates background from 

130 A® — > (A — > p7r~)// + /i~ which peaks above the B mass. For the B° — > K*°n + jj,^ decay, the 

131 same peaking background vetoes are used as in Ref. [15J, which remove contaminations 

132 from B° s — > <pfi + n~ , B° — > J/ip K*° and B° — > K*°fi + fi~ decays where the kaon and pion 

133 are swapped. Finally, for the B + — > K + fi + fi~ decay, backgrounds from B + — > J/ipK + 

134 and B — > ip(2S)K + are present, where the K + and fi + candidates are swapped. If a 

135 candidate has a K + fi~ track combination consistent with originating from a J/ip or 

136 ip(2S) resonance, the kaon is required to be inside the acceptance of the muon system 

137 but to have insufficient hits in the muon stations to be classified as a muon. These vetoes 

138 remove less than 1% of the signal and reduce peaking backgrounds to a negligible level. 

139 

wo The mass distribution of B candidates is shown versus the di-muon mass for B + — > 

141 K + fi + fi~ data in Fig. [TJ The other signal channels have similar distributions, but with 

142 a smaller number of events. The excess of candidates seen as horizontal bands around 

143 3090MeV/c 2 and 3690MeV/c 2 are due to J/ip and ip{2S) decays, respectively. These 

144 events are removed from the signal channels by excluding the di-muon regions in the 

145 ranges 2946— 3181 MeV/c 2 and 3586— 3766 MeV/c 2 . If a B candidate has an mass below 

146 5220MeV/c 2 the veto is extended to 2800-3181 MeV/c 2 and 3450-3766 MeV/c 2 to elim- 

147 inate candidates for which the J/ip or the ip(2S) decay undergoes final state radiation. 

148 Such events are shown in Fig. [I]as regions (a). In a small fraction of events, the di-muon 

149 mass is poorly reconstructed. This causes the J/ip and ip{2S) decay to leak into the re- 

150 gion just above the B mass. These events are shown in Fig. [l] as regions (b). The veto 

151 is extended to 2946— 3250 MeV/c 2 and 3586— 3816 MeV/c 2 in the candidate B mass region 

152 from 5330— 5460 MeV/c 2 to eliminate these events. These vetoes largely remove the char- 

153 monium resonances and reduce the combinatorial background. Regions (c) in Fig. [T] are 

154 composed of B — > J/ipK + X and B — > ip(2S)K + X decays where X is not reconstructed. 

155 In the subsequent analysis only candidates with masses above 5170MeV/c 2 are included 
we to avoid dependence on the shape of this background. 
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Figure 1: Mass of the di-muon versus the mass of the B + — > K + fi + fi~ candidates. Only 
the di-muon mass region close to the J/ip and ip(2S) masses is shown. The lines show the 
boundaries of the regions which are removed. Regions (a)-(c) are explained in the text. 

Table 1: Signal yields of the B— > K^*' fi + fi~ decays. The upper bound of the highest q 2 bin, 
g^ ax , is 19.3GeV 2 /c 4 and 23.0GeV 2 /c 4 for B-^ K*fx+^~ and B^ Kfj,+^~, respectively. 



q 2 range 




V 


fsT + /i + /i- 








[GeV 2 /c 4 ] 


L 


D 


L + D 


L 


D 


L + D 


0.05- 2.00 


1 ±2 


2±3 


135 ±13 


4±3 


5±4 


108 ± 11 


2.00- 4.30 


2±3 


-1 ±3 


175 ±16 


3 ±2 


5±3 


53 ± 9 


4.30- 8.68 


9±4 


16 ±6 


303 ± 22 


4±3 


17 ±6 


203 ± 17 


10.09-12.86 


4±3 


10 ±4 


214 ±18 


4±3 


15 ±5 


128 ± 14 


14.18 - 16.00 


3±2 


3±3 


166 ±15 


5±3 


4±3 


90 ± 10 


16.00- q 2 meix 


5±3 


4±3 


257 ± 19 


2 ± 1 


4±3 


80 ± 11 


1.00- 6.00 


8±4 


3±6 


356 ± 23 


5±3 


15 ±5 


155 ± 15 


0-05- q 2 max 


25 ±8 


35 ±11 


1250 ± 42 


23 ±6 


53 ± 10 


673 ± 30 



is? 4 Signal yield determination 

158 The yields for the signal channels are determined using extended unbinned maximum like- 

159 lihood fits to the ji + mass in the range 5170-5700 MeV/c 2 . These fits are performed 
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Figure 2: Mass distributions and fits of the signal channels integrated over the full q 2 
region. For the K® channels, the plots are shown separately for the L and D re- 
construction categories, (a,b) and (c,d) respectively. The signal component is shown by 
the dashed line, the partially reconstructed component in 2(a) and 2(c) is shown by the 
dotted line while the solid line shows the entire fit model. 



wo in seven q 2 bins and over the full range as shown in Table [l| The results of the fits inte- 

161 grated over the full q 2 range are shown in Fig. [2] After selection, the mass of K® candidates 

162 is constrained to the nominal K® mass. The signal component is described by the sum 
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163 of two Crystal Ball functions [TU] with common peak and tail parameters, but different 

164 widths. The shape is taken to be the same as the B— > J/ipK^*' normalisation channels. 

165 The combinatorial background is fitted with a single exponential function. As stated 
lee in Sect. |3j part of the combinatorial background is removed by the charmonium vetoes. 

167 This is accounted for by scaling the remaining background. For the B — >• Kfi + [i~ decays, 

168 a component arising mainly from partially reconstructed B— > K* / u + decays is present 
leg at masses below the B mass. This partially reconstructed background is characterised 

170 using a threshold model detailed in Ref. [IT]. The shape of the partial reconstruction 

171 component is again assumed to be the same as for the normalisation channels. For the 

172 B + — > K + fi + fi" channel, the impact of this component is negligible due to the relatively 

173 high signal and low background yields. For the B° — > Kgfi + [i~ channel, the amount of 

174 partially reconstructed decays is found to be less than 25% of the total combinatorial 

175 background in the fit range. 

176 The signal-shape parameters are allowed to vary in the B° — > J/ip K® mass fits 

177 and are subsequently fixed for the B° — > K®fi + fi~ mass fits when calculating the 

178 significance. The significance a of a signal S for B° — > K®fi + fi~ is defined as 

179 a 2 = 21n£ L (S) + 21n£ D (S)-21n£ L (0)-21n£ D (0) where C L > D (S) and £ L < D (0) are the 
lso likelihoods of the fit with and without the signal component, respectively. The 

181 B°—> Kgfj, + n~ channel is observed with a significance of 5.7 a. 

182 5 Normalisation 

183 In order to simplify the calculation of systematic uncertainties, each signal mode is nor- 

184 malised to the B — > J/ijj channel, where the J/tp decays into two muons. These decays 

185 have well measured branching fractions which are approximately two orders of magnitude 
lee higher than those of the signal decays. Each normalisation channel has similar kinematics 
187 and the same final state particles as the signal modes. 

las The relative efficiency between signal and normalisation channels is estimated using 

189 simulated events. After smearing the IP resolution of all tracks by 20%, the IP distribu- 

190 tions of candidates in the simulation and data agree well. The performance of the PID is 

191 studied using the decay D* + — > (D°— > 7t + K^)tt + , which provides a clean source of kaons 

192 to study the kaon PID efficiency, and a tag-and-probe sample of B + — > J/ipK + to study 

193 the muon PID efficiency The simulation is reweighted to match the PID performance of 

194 the data. 

195 Integrating over q 2 , the relative efficiency between the signal and normalisation chan- 
ge nels is between 70 and 80% depending on the decay mode and category. The relative 

197 efficiency includes differences in the geometrical acceptance, as well as the reconstruction, 

198 selection and trigger efficiencies. Most of these effects cancel in the efficiency ratio be- 

199 tween K® and K + channels, as shown in Fig. |3} The dominant effect remaining is due 

200 to the Kg reconstruction efficiency, which depends on the momentum. At low q 2 , the 

201 efficiency for B° — > K®[i + [i~ (D) decreases with respect to that for B + — > K + ji + \i~ due to 

202 the high K® momentum in this region. This results in the K® meson more often decaying 
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Figure 3: Efficiency of the K® channels with respect to the K + channels for (left) 
B^-Kfi + fi~ and (right) B^-K*fi + fi~, calculated using the simulation. The efficien- 
cies are shown for both L and D K® reconstruction categories and include the visible 
branching fraction of K° — » K® — > tt + tt~. The error bars are not visible as they are 
smaller than the marker size. 

203 beyond the TT and consequently it has a lower reconstruction efficiency This effect is 

204 not seen in the B + — > K* + fi + fi~ D category as the K® typically has lower momentum in 

205 this decay and so the K® reconstruction efficiency is approximately constant across q 2 . 

206 This Kg reconstruction effect is also seen in the L category for both modes but is par- 

207 tially compensated by the fact that the K® daughters can cause the event to be triggered, 

208 which increases the trigger efficiency with respect to the K + channels at low q 2 . Summed 

209 over both the L and D categories, the efficiency of the decays involving a K° meson is 

210 approximately 10% with respect to those involving a charged kaon. This is partly due to 
2n the visible branching fraction of K° — > K® — > 7t + ?t~ (~30%) and partly due to the lower 

212 reconstruction efficiency of the K® due to the long lifetime and the need to reconstruct an 

213 additional track (~30%). The relative efficiency between the L and D signal categories is 

214 cross-checked by comparing the ratio for the B ip(2S)K^ decay to the corresponding 

215 ratio for the B — > J/ip K^*> decays seen in data. The results agree within the statistical 

216 accuracy of 5%. 

217 6 Systematic uncertainties 

218 Gaussian constraints are used to include all systematic uncertainties in the fits for A\ 

219 and the branching fractions. In most cases the dominant systematic uncertainty is that 

220 from the branching fraction measurements of the normalisation channels, ranging from 3 

221 to 6%. There is also a statistical uncertainty on the yield of the normalisation channels, 

222 which is in the range 0.5-2.0%, depending on the channel. 

223 The finite size of the simulation samples introduces a statistical uncertainty on the 

224 relative efficiency and leads to a systematic uncertainty in the range 0.8-2.5% depending 
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225 on q 2 and decay mode. 

226 The relative tracking efficiency between the signal and normalisation channels is cor- 

227 rected using data. The statistical precision of these corrections leads to a systematic 

228 uncertainty of ~ 0.2% per long track. The differences between the downstream tracking 

229 efficiency between the simulation and data are expected to mostly cancel in the normal- 

230 isation procedure. A conservative systematic uncertainty of 1% per downstream track is 

231 assigned for the variation across q 2 . 

232 The PID efficiency is derived from data, and its corresponding systematic uncertainty 

233 arises from the statistical error associated with the PID efficiency measurements. The 

234 uncertainty on the relative efficiency is determined by randomly varying PID efficiencies 

235 within their uncertainties, and recomputing the relative efficiency. The resulting uncer- 

236 tainty is found to be negligible. 

237 The trigger efficiency is calculated using the simulation. Its uncertainty consists of 

238 two components, one associated with the trigger efficiency of the K® meson, and one 

239 associated with the trigger efficiency of the muons (and pion from the K*). For the 

240 muons and pion the uncertainty is obtained using B + — > J/ipK + and B° — > J/ipK*° 

241 events in data that are triggered independently of the signal. These candidates are used 

242 to calculate the trigger efficiency and are compared to the efficiency calculated using the 

243 same method in simulation. The difference is found to be ~ 2% for both B + — > J/tp K + 

244 and B° — > J/ip K*° decays and is assigned as a systematic uncertainty. This uncertainty is 

245 assumed to cancel for the isospin asymmetry as the presence of muons is common between 

246 the K® channels and the K + channels. The uncertainty associated with the K® trigger 

247 efficiency is calculated by comparing the fraction of candidates triggered by daughters 

248 in the simulation and the data. The difference is used as an estimate of the capability of 

249 simulation to reproduce these trigger decisions. The simulation is found to underestimate 

250 the Kg trigger decisions by 10-20% depending on the decay mode. This percentage is 

251 multiplied by the fraction of trigger decisions where the K® participates in a given bin of 

252 q 2 leading to an uncertainty of 0.2-4.1% depending on q 2 and decay mode. 

253 The effect of the unknown angular distribution of B + — > K* + fi + {i~ decays on the 

254 relative efficiency is estimated by altering the Wilson coefficients appearing in the operator 

255 product expansion method [TH ITU] . The Wilson coefficients, C7 and C10, have their real 

256 part inverted and the relative efficiency is recalculated. This can be seen as an extreme 

257 variation which is used to obtain a conservative estimate of the associated uncertainty. 

258 The calculation was performed using an EvtGen physics model which uses the transition 

259 form factors detailed in Ref. [20]. The difference in the relative efficiency varies from 0-6%, 

260 depending on q 2 , and it is assigned as a systematic uncertainty. 

261 The shape parameters for the signal modes are assumed to be the same as the nor- 

262 malisation channels. This assumption is validated using the simulation and no systematic 

263 uncertainty is assigned. The statistical uncertainties of these shape parameters are prop- 

264 agated through the fit using Gaussian constraints, accounting for correlations between 

265 the parameters. The uncertainty on the amount of partially reconstructed background 

266 is also added to the fit using Gaussian constraints, therefore no further uncertainty is 

267 added. The parametrisation of the fit model is cross-checked by varying the fit range 
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268 and background model. Consistent yields are observed and no systematic uncertainty is 

269 assigned. 

270 Overall the systematic error on the branching fraction is 4-8% depending on q 2 and 

271 the decay mode. This is small compared to the typical statistical error of ~ 40%. 

272 7 Results and conclusions 

273 The differential branching fraction in the i th q 2 bin can be written as 

d& _ N\B-+ KMfi+iT) B(B-> J/ip KW)B(J/iIj ->■ /i+aO 

~df~ N{B^J/^K(*)) X ' ^' 

274 where N l (B^ HfW 

fi + fi ) is the number of signal candidates in bin i, 

275 is the number of normalisation candidates, the product of B(B—y J/ipK^) and 

276 B(J/tp — > is the visible branching fraction of the normalisation channel [21], e* el is 

277 the relative efficiency between the signal and normalisation channels in bin i and finally 

278 A* is the bin i width. The differential branching fraction is determined by simultane- 

279 ously fitting the L and D categories of the signal channels. The branching fraction of 

280 the signal channel is introduced as a fit parameter by re-arranging Eq. (|2]) in terms of 

281 N(B—y K^*> fi + n~). Confidence intervals are evaluated by scanning the profile likelihood. 

282 The results of these fits for B° — y K°fi + pT and B + — y K* + fi + fx~ decays are shown in 

283 Fig. [4] and given in Tables [2] and [3j Theoretical predictions [221 1213 [21] are superimposed 

284 on Figs. [4] and [HJ In the low q 2 region, these predictions rely on the QCD factorisation 

285 approaches from Refs. [251 [22] for 5— y K*ji + n~ and Ref. [27J for B—y Kn + n~ which lose 

286 accuracy when approaching the J/0 resonance. In the high q 2 region, an operator product 

287 expansion in the inverse 6-quark mass, l/m&, and in \j \fq 2 is used based on Ref. [28] . 

288 This expansion is only valid above the open charm threshold. In both q 2 regions the 

289 form factor calculations for B — y K*fi + fi~ and B —y Kfi + fi~ are taken from Refs. [2H] 

290 and [30J respectively. These form factors lead to a high correlation in the uncertainty 

291 of the predictions across q 2 . A dimensional estimate is made of the uncertainty from 

292 expansion corrections [31]. The non-zero isospin asymmetry arises in the low q 2 region 

293 due to spectator-quark differences in the so-called hard-scattering part. There are also 

294 sub-leading corrections included from Refs. [1] and |26j which only affect the charged 

295 modes and further contribute to the isospin asymmetry. 

296 The total branching fractions are also measured by extrapolating underneath the char- 

297 monium resonances assuming the same q 2 distribution as in the simulation. The branching 

298 fractions of B° — y K°/i + fi~ and B + —y K* + fi + fi~ are found to be 

B(B ^ K°fi + fi-) = (0.3lj£{£) x 10~ 6 and 
B(B + ^ K* + fi + i2~) = (1.16 ± 0.19) x 10~ 6 , 

299 respectively, where the errors include statistical and systematic uncertainties. These re- 

300 suits are in agreement with previous measurements and with better precision |21j . 
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Figure 4: Differential branching fractions of (left) B° — > K°fi + fi and (right) 
B + — > K* + fi + fi~. The theoretical SM predictions are taken from Refs. 
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Figure 5: Isospin asymmetry of (left) B^-Kfi + fi and (right) B^-K*fi + ^ . For 
B— > K*fi + fi~ the theoretical SM prediction, which is very close to zero, is shown for 
q 2 below 8.68GeV/c 2 , from Ref. [21]. 



301 The isospin asymmetries as a function of q 2 for B — > K[x + fi~ and B — > K*fi + fi~ are 

302 shown in Fig. [5] and given in Tables [2] and [3j As for the branching fractions, the fit is done 

303 simultaneously for both the L and D categories where Aj is a common parameter for the 

304 two cases. The confidence intervals are also determined by scanning the profile likelihood. 

305 The significance of the deviation from the null hypothesis is obtained by fixing A\ to be 

306 zero and computing the difference in the negative log-likelihood from the nominal fit. 

307 In summary, the isospin asymmetries of B — y K^fj, + fx~ decays and the branching 

308 fractions of B° — > K°fi + fi~ and B + — > K* + fi + fi~ are measured, using 1.0 fb _1 of data taken 

309 with the LHCb detector. The two q 2 bins below 4.3GeV/c 2 and the highest bin above 

310 16 GeV/c 2 have the most negative isospin asymmetry in the B—y Kfi + fi~ channel. These 

311 q 2 regions are furthest from the charmonium regions and are therefore cleanly predicted 
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Table 2: Partial branching fractions of B°^-K°jj + j2~ and isospin asymmetries of 
B— > K[i + pT decays. The significance of the deviation of A\ from zero is shown in the 
last column. The errors include the statistical and systematic uncertainties. 



q 2 range [ GeV 2 /c 4 ] dB/dq 2 [10" 8 / GeV 2 / c 41 



A, 



a(Ax = 0) 



0.05- 2.00 
2.00- 4.30 
4.30- 8.68 
10.09-12.86 
14.18 - 16.00 
16.00 - 23.00 



1 1 +L4 

L - L -1.2 

3 +L1 
2.8 ±0.7 

1 O+0.8 
i -°-0.7 

1 l +0 - 7 

5+°- 3 

U -°-0.2 



-0.55^; 4 ° 1.5 



-0.76 



H).45 
-0.79 



o.oor^ 4 



-0.40 ±0.22 



1.9 
0.1 



-0.15±g;il 0.8 



1.9 



-0.52^ 3.0 



1.00- 6.00 



i q+0.9 
1 -°-0.7 



-0.35±g;|? 1.7 



Table 3: Partial branching fractions of B + — > K* + fi + n~ and isospin asymmetries of 
B— > K*fi + fi~ decays. The significance of the deviation of A\ from zero is shown in the 
last column. The errors include the statistical and systematic uncertainties. 



q 2 range [ GeV 2 /c 4 ] dB/dq 2 [10" 8 / GeV 2 /c 4 



-4, 



a(A 1 = 0) 



0.05- 2.00 
2.00- 4.30 
4.30- 8.68 
10.09-12.86 
14.18 - 16.00 
16.00-19.30 



5.7tl- 7 

7 7+2.6 
'■ '-2.4 

5 5+ 2 ' 6 
3.8 ± 1.4 



-o.o6i$;i5 



0.02^; 2 ? 
02+ - 21 

u - uz -0.20 



0.2 



-0.27t° il 0.9 



0.4 



-0.16l^ 6 7 0.9 



0.1 
0.1 



1.00- 6.00 



5.81J; 8 



-0.15 ±0.16 



1.0 



312 theoretically. This asymmetry is dominated by a deficit in the observed B° — > K°/j, + jj,^ 

313 signal. Ignoring the small correlation of errors between each q 2 bin, the significance of 

314 the deviation from zero integrated across q 2 is calculated to be 4.4 o. The B— > K*fi + fi~ 

315 case agrees with the SM prediction of almost zero isospin asymmetry [I]. All results agree 

316 with previous measurements [31 [321 [33] . 



317 Acknowledgements 

318 We would like to thank Christoph Bobeth, Danny van Dyk and Gudrun Hiller for provid- 

319 ing SM predictions for the branching fractions and the isospin asymmetry of B — > K*fi + fi~ 

320 decays. We express our gratitude to our colleagues in the CERN accelerator departments 

321 for the excellent performance of the LHC We thank the technical and administrative 



12 



322 staff at CERN and at the LHCb institutes, and acknowledge support from the Na- 

323 tional Agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); CERN; NSFC (China); 

324 CNRS/IN2P3 (France); BMBF, DFG, HGF and MPG (Germany); SFI (Ireland); INFN 

32 5 (Italy); FOM and NWO (The Netherlands); SCSR (Poland); ANCS (Romania); MinES 

326 of Russia and Rosatom (Russia); MICINN, XuntaGal and GENCAT (Spain); SNSF and 

327 SER (Switzerland); NAS Ukraine (Ukraine); STFC (United Kingdom); NSF (USA). We 

328 also acknowledge the support received from the ERC under FP7 and the Region Au- 

329 vergne. 



330 References 

331 [1] T. Feldmann and J. Matias, Forward-backward and isospin asymmetry for B — > 

332 K(*H + l~ decay in the Standard Model and in super symmetry, JHEP 01 (2003) 074, 

333 |arXiv : hep-ph/02121581 

334 [2] BaBar collaboration, B. Aubert et al., Direc t CP, l epton flavor an d isospin asym- 

335 metries in the decays B K^l + l~ , |Phys. Rev. Lett. 102 (2008) 0918031 

336 larXiv:0807.41191 

337 [3] Belle collaboration, J.-T. Wei et al., Measurement of the differential branching frac- 

338 tion and forward-backward asymmetry for B — > K^*H + l~ , \Phys. Rev. Lett. 103 (200*9)] 

339 |171801j [arXiv : 0904 . 07701 



340 



[4] Heavy Flavor Averaging Group, D. Asner et al., Averages of b-hadron, c-hadron, and 
34i r -lepton properties, arXi v: 1010. 1589. 



342 [5] LHCb collaboration, A. A. Alves Jr. et al, The LHCb detector at the LHC, JINST 

343 |3 (2008) S 08005 



344 n [6] V. Gligorov, C. Thomas, and M. Williams, The HLT inclusive B triggers, LHCb- 

345 IPUB-20TT- 016, 

346 [7] T. Sjostrand, S. Mrenna, and P. Skands, PYTHIA 6.4 Physics and manual, [JHEP 

347 |05 (2006) 02"6l|arXiv:hep-ph/0603175l 

348 [8] I. Belyaev et al, Handling of the generation of primary events in GAUSS, the LHCb 



349 simulation framework, Nuclear Science Symposium Conference Record (NSS/MIC) 

|IEEE (2010) 1155| 



351 _ [9] D. J. Lange, The EvtGen particle decay simulation package, Nucl. Instrum. Meth. 



352 



A462 (2001) 152. 



353 [10] P. Golonka and Z. Was, PHOTOS Monte Carlo: a precision tool for QED corrections 

354 in Z and W decays, |Eur. Phys. J. C45 (2006) 97| |arXiv : hep-ph/0 506026 , 



13 



355 [11] A. Ali, P. Ball, L. T. Handoko, and G. Hiller, Comparative study of the decays B — > 

356 (K , K*)£ + £~ in Standard Model and super symmetric theories, |Phys. Rev. D61 ( 2000) 

357 10740241 [arXiv : hep-ph/9910221[ 

358 [12] GEANT4 collaboration, J. Allison et al., GeantJ^ developments and applications, 

359 |IEEE Trans. Nucl. Sci. 53 (2006) 270] GEANT4 collaboration, S. Agostinelli et al, 



see GEANT4: A simulation toolkit, Nucl. Instrum. Meth. A506 (2003) 250 



361 [13] M. Clemencic et al, The LHCb simulation application, Gauss: design, evolution and 

362 experience, |J. of Phys: Conf. Ser 331 (2011) 032023] 

363 [14] P. Speckmayer, A. Hocker, J. Stelzer, and H. Voss, The toolkit for multivariate data 

364 analysis TMVA 4, |X~oTPhys: Conf. Ser 219 (2010) 032057| 

365 [15] R. Aaij et al, Differential branching fraction and angular analysis of the decay B° — > 

366 K* i2 + n~, Phys. Rev. Lett. 108 (2012) 181806, |arXiv: 1112. 3515} 

367 [16] T. Skwarnicki, A study of the radiative cascade transitions between the Upsilon-prime 

368 and Upsilon resonances, PhD thesis, Institute of Nuclear Physics, Krakow, 1986, 
IDESY-F31-8F021 



370 [17] P. Koppenburg, Selection of B + — > UK at LHCb and sensitivity to Rk, LHCb-2007- 

371 IU341 



372 [18] K. G. Wilson, Non-Lagrangian models of current algebra, Phys. Rev. 179 (1969) 

373 H1991 

374 [19] K. G. Wilson and W. Zimmermann, Operator product expansions and composite field 

375 operators in the general framework of quantum field theory, Commun. Math. Phys. 

376 |24 (1972) 87] 



377 [20] D. Melikhov, B. Stech, Weak form factors for heavy meson decays: an update, Phys 



378 



379 



Rev. D62 (2000) 014006, arXiv:hep-ph/0001113 



[21] Particle Data Group, K. Nakamura et al, Review of particle physics, J. Phys. G37 

380 1(2010) 07502T| 

381 [22] C. Bobeth, G. Hiller, and D. van Dyk, More benefits of semileptonic rare B decays 



382 at low recoil: CP Violation, JHEP 07 (2011) 067, arXiv : 1105 . 0376 



383 [23] C. Bobeth, G. Hiller, D. van Dyk, and C. Wacker, The Decay B K£ + £~ at low 



384 hadronic recoil and model-independent AB = 1 constraints, JHEP 01 (2012) 107 



385 arX iv: 1111.2 558, 

386 [24] Private communication with C. Bobeth, G. Hiller and D. van Dyk. 

387 [25] M. Beneke, T. Feldmann, and D. Seidel, Systematic approach to exclusive B — > 

388 Vl+l-, V-f decays, [Nucl. Phys. B612 (2001) 25] |arXiv:hep-ph/0106067| 



14 



M. Beneke, T. Feldmann, and D. Seidel, Exclusive radiative and electroweak 
b — >• d and b — >• s penguin decays at NLO, |Eur. Phys. J. C41 (2005) 173[ 
|arXiv : hep-ph/04 12400} 

C. Bobeth, G. Hiller, and G. Piranishvili, Angular distributions of B — >■ Kit decays, 
|JHEP 12 (2007) 040|[^-Xiv: 0709 .41741 

B. Grinstein and D. Pirjol, Exclusive rare B — > K*£ + £~ - decays at low re- 
coil: controlling the long-distance effects, |Phys. Rev. D70 (2004) 114 005, 



arXiv : hep-ph/0404250 , 

P. Ball and R. Zwicky, B^^ — > P, K* , <fi decay form factors from light-cone sum 
rules reexamined, |Phys. Rev. D71 (2005) 014029[ |arXiv : hep-ph/04 12079) 

A. Khodjamirian, T. Mannel, A. Pivovarov, and Y.-M. Wang, Charm-loop effect in 
B K^i + i' and B -» K*j, |JHEP 09 (2010) 089[ |arXiv: 1006 . 4945| 



U. Egede et al, New observables in the decay mode B° d -»■ K*°£ + £~, |jHEP 11 (2008)] 
[0H2l|arXiv: 0807 .2589| 

BaBar collaboration, Measurement of branching fractions and rate asymmetries in 



the rare decays B -> K^£ + £~, submitted to Phys. Rev. D, arXiv: 1204.3933 



CDF collaboration, T. Aaltonen et al., Observation of the baryonic flavor- changing 



neutral current decay A° b A[i + [i~ , |Phys. Rev. Lett. 107 (2011) 20 1802 

larXiv: 1107.37531 



15 



